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Abstract 

A molecular model for carbon dioxide is assessed regarding vapor-liquid equilibrium properties. 
Large deviations, being above 15 %, are found for vapor pressure and saturated vapor density in 
the entire temperature range. 
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In a recent article, Zhang and Duan^ presented a new potential model for carbon dioxide 
(CO2). It consists of three Lennard- Jones (LJ) sites to account for repulsion and dispersion 
and three distributed partial charges to describe the quadrupolar interaction. The molecular 
model is rigid and rotationally symmetric around the molecular axis. In that work^, simula- 
tion results on vapor-liquid equilibria (VLE), radial distribution function and self-diffusion 
coefficient have been reported for the new CO2 model which are in excellent agreement with 
experimental data. They also compared to results from different other models from the 
literature. The new model is found "to be superior to the previous models in general". 

For the VLE properties, deviations between model and experimental data are reported to 
bei: 0.7 % for vapor pressure, 0.1 % for saturated liquid density, 2.3 % for saturated vapor 
density, and 1.9 % for heat of vaporization. 

Particularly for vapor pressure, which is the most sensitive of those properties, the authors 
claim to achieve an average accuracy of better than 1 % over the entire temperature range 
with a molecular model which is noteworthy for any molecular model. 

As there is widespread scientific interest in CO2, the model by Zhang and Duan^ was 
employed for subsequent work-'^'^'^i^. It should be pointed out that no VLE data were 
published in^i^iii^i^. We recently tested the CO2 model from Zhang and Duan^ and found 
results that strongly deviate from these reported by Zhang and Duani. The new simulation 
results also strongly deviate from experimental data, particularly the vapor pressure and 
the saturated vapor density. 

The present assessment was made on the basis of two simulation tools that employ dif- 
ferent methods to determine VLE. Firstly, the Grand Equilibrium method^ was used as 
implemented in our simulation tool ms2^ and, secondly, the Gibbs ensemble^ was used as 
implemented in the freely available simulation tool TOWHEE^^. Both programs have proven 
to be correct e.g. in the recent Industrial Fluid Property Simulation Challenge for the case 
of ethylene oxide which is very similar to the present case^. 

For our approach with ms2, molecular dynamics (MD) simulations were performed in 
the liquid phase containing 1024 molecules. After a sufficient equilibration period, the 
chemical potential was calculated by Widoms's insertion method^^ over 300 000 time steps. 
According to the Grand Equilibrium method^,^ the dew point was sampled with Monte-Carlo 
(MC) simulations where approximately 500 molecules were used. The simulation details were 
similar to those published and are not repeated here. 
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With TOWHEE^°, Gibbs ensemble MC simulationa^ were performed. There, smaller 
systems were studied, containing 500 molecules. After an equilibration over 5 000 loops 
without volume and molecules transfer moves, followed by 10 000 loops with these moves, 
50 000 production loops were performed. Other simulation details were similar to those 
published in^^ and are not repeated here. 

The Present simulation results are compared to those reported by Zhang and Duan^ 
in Table [T] as well as in Figures [1] and [2|, which also contains results from a highly accurate 
reference data for CO-r^ recommended by the National Institute of Science and Technology^^. 
Figure 2 shows deviation plots where it can be seen that the present data sets based on the 
two different simulation methods agree with each other within their (combined) error bars 
throughout. However, they are significantly off the data by Zhang and Duan^ which coincide 
excellently with the experiment. 

The present simulation data (Grand Equilibrium) show the following average deviations 
from experimental data: 18 % for vapor pressure, 0.6 % for saturated liquid density, 17 % 
for saturated vapor density and 4.6 % for heat of vaporization. Only the saturated liquid 
density by Zhang and Duan^ and the heat of vaporization are in good agreement with the 
experimental data, the vapor pressure and saturated vapor density are significantly too high 
throughout most of the temperature range. 

It has to be concluded that the CO2 model by Zhang and Duan^ is not generally superior 
to previous models. 
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Figure 1 

Vapor pressure of carbon dioxide. • present results Grand Equilibrium method'', A present 
results Gibbs ensemble method^, ■ results by Zhang and Duan^, — reference EOS^. 

Figure 2 

Relative deviations of vapor-liquid equilibrium properties between simulation and reference 
EOS^* {5z = {zsim — -2^eos)/-2^eos)- • prcscut rcsults Grand Equilibrium method"^, A present 
results Gibbs ensemble method-, ■ results by Zhang and Duan^, — reference EOS^. 
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TABLE I: Vapor-liquid equilibria of carbon dioxide: present simulation results with the Grand Equilibrium method (GE)^ and the Gibbs 
ensemble (Gibbs)^ compared to simulation results by Zhang and Duan (Zhang)^ and the reference EOS (eos)^. The number in parentheses 
indicates the statistical uncertainty in the last digit. 



T 


PGE 


PGibbs PZhang 


Pcos 


p'ge 


^'Gibbs Pzhang 


Pcos 


p'6b 


^Gibbs ^Zhang 


p" 

rcos 




A/^Sibbs 


^"■Zhang 


^Kos 


K 


MPa 


MPa MPa 


MPA 


mol/1 


mol/1 mol/1 


mol/1 


mol/1 


mol/1 mol/1 


mol/1 


kJ /mol 


k J / mol 


kJ / mol kJ /mol 


220 


0.77(1) 




0.599 


26.42(3) 




26.497 


0.466(1) 




0.359 


14.41 (1) 






15.179 


230 


1.08(1) 


— 0.90 (2) 


0.893 


25.48(2) 


~ 25.59 (1) 


25.646 


0.641(2) 


— 0.558(7) 


0.529 


13.717(6) 




13.90(5) 


14.437 


240 


1.50(2) 


— 1.29 (1) 


1.283 


24.60(2) 


— 24.697(6) 


24.742 


0.882(2) 


— 0.79 (2) 


0.757 


12.994(7) 




13.19(6) 


13.628 


250 


2.11(2) 


— 1.788(2) 


1.785 


23.66(3) 


— 23.740(8) 


23.767 


1.253(4) 


— 1.09 (3) 


1.06 


12.138(8) 




12.40(6) 


12.733 


260 


2.86(2) 


2.9(2) 2.41 (2) 


2.419 


22.58(4) 


22.6(2) 22.68 (1) 


22.697 


1.738(6) 


1.8(1) 1.48 (4) 


1.464 


11.14 (1) 


11.2(1) 


11.51(5) 


11.728 


270 


3.69(2) 


3.8(3) 3.18 (3) 


3.203 


21.33(6) 


21.3(3) 21.50 (1) 


21.491 


2.29 (1) 


2.4(3) 2.01 (5) 


2.008 


10.06 (1) 


10.0(2) 


10.46(4) 


10.569 


280 


4.78(3) 


4.8(2) 4.11 (5) 


4.161 


19.9 (1) 


20.1(2) 20.09 (2) 


20.077 


3.14 (2) 


3.2(2) 2.76 (4) 


2.766 


8.71 (2) 


8.8(1) 


9.17(3) 


9.183 


290 


5.92(4) 


— 5.23 (6) 


5.318 


18.1 (4) 


— 18.29 (4) 


18.284 


4.14 (4) 


— 3.96 (6) 


3.907 


7.16 (4) 




7.45(4) 


7.399 



